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Non-Gaussian velocity shears 
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Abstract. We report on a novel kind of small scale structure in molecular clouds found in IRAM-30m and CSO 
maps of 12 CO and 13 CO lines around low mass starless dense cores. These structures come to light as the locus 
of the extrema of velocity shears in the maps, computed as the increments at small scale (~ 0.02 pc) of the line 
velocity centroids. These extrema populate the non-Gaussian wings of the shear probability distribution function 
(shear-PDF) built for each map. They form elongated structures of variable thickness, ranging from less than 0.02 
pc for those unresolved, up to 0.08 pc. They are essentially pure velocity structures. We propose that these small 
scale structures of velocity shear extrema trace the locations of enhanced dissipation in interstellar turbulence. In 
this picture, we find that a significant fraction of the turbulent energy present in the field would be dissipating in 
structures filling less than a few % of the cloud volume. 
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1. Introduction 

Star formation proceeds at vastly different rates, in space 
and time, within a given galaxy and from one galaxy to 
another. These rates are known to be governed by the lo- 
cal conditions prevailing in dense and cold gas, but also 
depend on large scale environments, up to extragalactic 
scales. The only two processes, together with rotation, 
able to mitigate the effects of gravity are magnetic fields 
and supersonic turbulence because they involve energies of 
the same order of magnitude as the gravitational energy in 
the densest phases of the interstellar medium (ISM). More 
precisely, turbulent energy, because of its steep power 
spectrum, can stabilize the largest masses, first prone to 
gravitational instability, a property not shared by thermal 
energy which is scale-free (Panis & Perault 1998). On the 
other hand, it has also been proposed that supersonic tur- 
bulence might trigger star formation in shocks (Klessen 
et al. 2000). The possible role of turbulence in the star 
formation process is the motivation behind the plethora 
of recent studies on interstellar turbulence. 

Important insights to the field have been provided by 
the determination of the multi-scale properties of inter- 
stellar clouds and their comparison to the scaling laws 
of turbulence. A broad variety of statistical tools have 
been used: wavelet analysis (Gill & Henriksen 1990; 



Langer et al. 1993), A-variance (Bensch et al. 2001), 
auto-correlation function (Kleiner & Dickman 1984, 1985; 
Dickman & Kleiner 1985; Perault et al. 1986), structure 
functions (Miesch & Bally 1994; Miesch & Scab 1995; 
Miesch et al. 1999; Padoan et al. 2003a), analysis in terms 
of fractal structures (Bazcll & Desert 1988; Falgarone 
et al. 1991; Stutzki et al. 1998). The principal compo- 
nent analysis has been used to diagnose the large-scale 
flows of atomic gas into which turbulence in molecular 
clouds is embedded (Brunt 2003). Observations of dense 
cores and their environment have revealed a break of the 
scaling properties of molecular clouds at the scale of the 
dense cores (Falgarone et al. 1998; Goodman et al. 1998). 
Direct numerical simulations of compressible turbulence 
have also been used to compare the real observables of the 
interstellar clouds to those simulated. The first attempt 
by Falgarone et al. (1994), using the high resolution sim- 
ulations of midly compressible turbulence of Porter et al. 
(1994) has been followed by more sophisticated compar- 
isons such as those based on the spectral correlation func- 
tion method (Rosolowsky et al. 1999). Numerical simula- 
tions of turbulence including magnetic fields, e.g. Ostriker 
et al. (2001), and self-gravity (Klessen et al. 2000; Heitsch 
et al. 2001) provided further support to the fact that in- 
terstellar turbulence bears many of the statistical proper- 
ties of supersonic magneto-hydrodynamical (MHD) tur- 
bulence, as simulated. 
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One major surprise brought to the field by direct nu- 
merical simulations of MHD turbulence has been that 
magnetic fields do not delay the dissipation of supersonic 
turbulence (Mac Low et al. 1998). This unexpected result 
calls for further numerical and observational approaches. 
The present paper is observational: It is an attempt to 
disclose kinematic signatures of turbulent dissipation in 
molecular clouds. 

To discuss what these kinematic signatures may be, we 
first need to briefly recall the main drivers of dissipation in 
interstellar turbulence, assuming infinite conductivity and 
thus neglecting Ohmic dissipation. The primary source of 
dissipation of supersonic turbulence is shocks, but shock 
interaction generates vorticity, as do non-planar shocks 
(Porter et al. 1994). Since the divergence of the veloc- 
ity field eludes direct detection in space, shocks cannot 
be detected by this kinematic signature, and remanent 
vorticity is a plausible signature of fossil shocks. Viscous 
dissipation is also present and is due primarily to elas- 
tic collisions. Dissipation due to neutral-neutral collisions 
follows the shear of the velocity field, and therefore the 
vorticity (Landau & Lifshitz 1987), and that due to ion- 
neutral collisions increases with the drift velocity of the 
ions relative to the neutrals (Kulsrud & Pearce 1969). 
This drift causes a force on the ions which, over timescales 
longer than the ion- neutral collision time (Zweibel 1988), 
is balanced by the Lorentz force (V x B) x B. The dissi- 
pation driven by ion-neutral collisions therefore involves 
J = V x B, the current density. 

Now, both vorticity, in hydrodynamical turbulence, 
and current density, in plasma turbulence, are known to 
be intermittent in space and time. Laboratory experiments 
of incompressible turbulence show the formation of tran- 
sient long and thin coherent vortices at the edge of which 
the velocity shear is so large that a significant fraction of 
the viscous dissipation occurs there (Douady et al. 1991). 
In plasma turbulence, intermittency also exists and is at 
the origin of non-Gaussian probability distribution func- 
tions (PDFs) of current density. It is seen in Tokamak 
plasma turbulence (Wang et al. 1999). Observations in 
the solar wind reveal an intermittent dissipation as well. 
The intermittency of dissipation of plasma turbulence has 
been modelled by Politano & Pouquet (1995), who pro- 
pose that the dissipative structures are sheetlike struc- 
tures of intense current density, in agreement with re- 
cent 3-dimensional numerical simulations of MHD turbu- 
lence (Politano et al. 1995; Biskamp & Miiller 2000). A 
recent attempt at estimating the sizescales of dissipation 
in incompressible MHD turbulence (Cho et al. 2002) shows 
that magnetic structures develop at scales much smaller 
than the viscous damping scale. 

Turbulence in molecular clouds is neither laboratory 
turbulence nor plasma turbulence, but the above elements 
suggest that dissipation follows the vorticity and the cur- 
rent density, and that these quantities are intermittent, 
i.e. exhibit large fluctuations at small scales. The origin of 
currents in molecular clouds is not known, but is likely as- 
sociated to differential rotation within clouds at the origin 



of toroidal or helical fields (Joulain et al. 1998). Such he- 
lical fields have been inferred from various observations of 
molecular clouds (Hanawa et al. 1993; Joulain et al. 1998; 
Carlqvist et al. 1998; Falgarone et al. 2001; Matthews et al. 
2002, 2001). These fields have been invoked to explain the 
polarization patterns of the dust continuum emission of 
filaments of matter (Harjunp et al. 1999; Fiege & Pudritz 
2000c) and their gravitational stability (Fiege & Pudritz 
2000a,b). Enhanced vorticity in molecular clouds should 
therefore be a good tracer of several major dissipative pro- 
cesses of turbulence in weakly ionized molecular clouds. Its 
kinematic signature, if available in spectral observations of 
molecular clouds, is essential to search for because it may 
trace fossil shocks, regions of large viscous dissipation in 
the neutrals, or local gas differential rotation at the origin 
of currents. 

In this paper, we analyze the velocity structure of sev- 
eral fields in molecular clouds in order to locate and char- 
acterize the regions of enhanced vorticity, likely to be re- 
gions of enhanced dissipation. We focus on the environ- 
ment of low mass dense cores almost thermally supported, 
where dissipation is anticipated to occur, or has occurred 
in a recent past. Section 2 gives the method used to trace 
vorticity. Section 3 is a description of the target fields. The 
statistical analysis of the velocity fields, and the manifesta- 
tions of their non-Gaussian features are given in Section 4. 
The spatial distribution of these non-Gaussian events is 
compared to that of the dense gas in Section 5. The pos- 
sible biases of the method are discussed in Section 6. The 
implications of this study in terms of turbulence dissipa- 
tion are discussed in Section 7. 

2. Method for measuring vorticity 

Measuring vorticity in interstellar turbulence is not 
straightforward because (1) the only velocity component 
accessible to measurement is its projection on the line of 
sight (v x ), and (2) measured spatial variations are lim- 
ited to those in the plane of the sky (y, z). We thus have 
access only to incomplete components of the vorticity 
(u>£ = dv x /d( where ( = y or z). By analyzing the outputs 
of numerical simulations of midly supersonic turbulence 
by Porter et al. (1994), Lis et al. (1996) have shown that 
vorticity extrema may be localized in a map of molecu- 
lar lines with high enough spectral resolution. Instead of 
computing the vorticity (an essentially inaccessible quan- 
tity in astronomy), they study statistical properties of line 
centroid velocity increments. More precisely, the centroid 
velocity (C) of a line is its first order moment 

C = J T(v x )v x dv x I J T(v x )dv x . (1) 

It is only in the optically thin case, that a line can be 
interpreted as the probability distribution function of the 
radial velocity along the observed line of sight. In this case, 
the centroid velocity is just the column density weighted 
average of the velocity along the line of sight. Centroid 
velocity increments associated to a lag I, are the differences 
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of centroid velocities for any pair of lines of sight separated 
by the distance I in the plane of sky. Thus, if r is the vector 
position and I is a vector of modulus I lying in the plane 
of sky, centroid velocity increments arc defined as 

5C{r,l) = C{r + l)-C{r). (2) 

Lis et al. (1996) show that the PDFs of centroid velocity 
increments have non-Gaussian wings at small lags, those 
wings gradually disappearing when the lag increases. More 
interestingly, they show that the spatial distribution of the 
positions populating these non-Gaussian wings is the same 
as that of the extrema of ((V x v) y ) 2 + ((V x v)^) 2 , where 
the brackets hold for the line-of-sight average. In the op- 
tically thin case, this result simply follows from the linear- 
ity of the operations of derivation, difference and summa- 
tion along a line of sight. The largest centroid increments 
thus trace a subset of the regions of largest vorticity, inte- 
grated along the line of sight. It has to be appreciated here 
that the averaging along the line of sight occurs first, over 
vorticity components which have a sign. An extremum of 
the above quantity means the occurrence along the line of 
sight of one (or a few) events far above the average. 

This method has been challenged by Klessen (2000) 
who computes the PDFs of line centroid increments (At>- 
PDFs) in numerical simulations of compressible turbu- 
lence with and without self-gravity. He finds that, when 
computed in simulations of decaying turbulence, Af -PDFs 
have shapes in disagreement with those observed. He 
shows that the inclusion of self-gravity in the simulations 
leads to better agreement with the Au-PDFs observed in 
molecular clouds. This discussion is valuable but its bear- 
ing is somewhat limited by the small size of the simula- 
tions (64 3 ). In addition, the confrontation with the obser- 
vations is restricted to star forming regions. It is therefore 
not surprising that the agreement with the numerical sim- 
ulations be better when self-gravity is included in the sim- 
ulations. Klessen concludes that "one should not rely on 
analyzing velocity PDFs alone to disentangle the different 
physical processes influencing interstellar turbulence" . We 
fully agree with him. Our scope here is far more limited. 
We search for regions of enhanced vorticity in turbulent 
interstellar clouds far from star forming regions, with the 
perspective of tracing bursts of dissipation of turbulence 
and characterizing them, on the basis of their statistical 
properties, spatial distribution and morphology. 

The spectral correlation function (SCF) introduced 
by Rosolowsky et al. (1999) is a general tool that is able to 
detect velocity and/or column density variations through 
an analysis of line shape differences. The flexibility of the 
SCF formalism makes it easy to use for direct compari- 
son between observations and simulations (Padoan et al. 
2003b). For instance, Ballesteros-Paredes et al. (2002) use 
the local form of the SCF to detect small scale veloc- 
ity structures in Hi observations of the north celestial 
pole loop. However, as recognized by those authors, this 
method cannot disentangle between the various origins 
of the small scale velocity variations {e.g. caustics versus 
shocks). The method we propose here is complementary as 



centroid velocity increments directly relate to the velocity 
shear (and thus vorticity) in the optically thin case. 

3. Description of the observed fields 

The fields analyzed in this paper fulfill several require- 
ments: (1) the observed maps are large enough to provide 
significant statistics on the velocity field, (2) the fields are 
far from star forming regions, in order not to include in- 
ternal energy sources (outflows, HII regions in expansion, 
...) in the statistics, (3) they comprise low mass, almost 
thermally supported dense cores. 

The two fields are nearby regions of low average 
column density, except in the small projection area of 
the dense cores. One, MCLD123.5+24.9 (hereafter re- 
ferred to as Polaris), is located in the high latitude cir- 
rus cloud of the Polaris Flare and the other, (hereafter 
L1512), is located at the eastern edge of the Taurus- Auriga 
complex. The average H2 column densities (at the arc 
min scale) deduced from 12 CO observations and/or star 
counts (Cambresy et al. 2001) in these fields are respec- 
tively N(R 2 ) ~ 10 21 and 2.5 xl0 20 cnr 2 or about 1 and 
0.25 mag of visual extinction, respectively. For compari- 
son, the visual extinction in the Taurus molecular cloud 
reaches 33 mag in the most opaque regions like TMC1 
or L1495 (Padoan et al. 2002), about one hundred times 
more opaque than the transparent areas that we analyze 
in this paper. 

These fields have been mapped in the two lowest ro- 
tational transitions of 12 CO and 13 CO as well as in C 18 
(J=l-0) at the IRAM-30m telescope (Falgarone et al. 
1998). A larger field around L1512 has been mapped 
with the 10.4 m antenna of the Caltcch Submillimeter 
Observatory (CSO) in the 12 CO (J=2-l) line and is de- 
scribed in Falgarone et al. (2001). The characteristics of 
the observations arc summarized in Table 1 and 2. Note 
that all the maps are Nyquist sampled or better. The aver- 
age signal-to-noise ratios (and their rms dispersion within 
the maps) are given in the last two columns of Table 2, for 
the line peak temperatures and the line integrated areas. 
These numbers illustrate the high quality of the data. The 
maps of the line integrated areas (zero order moment), and 
line centroid velocities (first order moment) are displayed 
in Figs. 1 through 3. We have inserted the high angu- 
lar resolution IRAM data at their relevant position in the 
large scale CSO field of L1512 to illustrate the agreement 
between the two data sets (Fig. 3). 

In spite of their low average extinction, each field har- 
bours a low mass dense core. The dense core in the Polaris 
field has been mapped in the CS(2-1), (3-2) and (5-4) lines 
(Hcithausen 1999) and in a number of molecular species at 
a few positions only by Gerin et al. (1997). Recent obser- 
vations of the dust thermal emission (Hcithausen et al. 
2002) show a clear peak of submillimeter dust contin- 
uum emission coinciding with the C 18 peak of Falgarone 
et al. (1998). Well-defined and barely resolved peaks of 
HC3N emission are found to be clearly displaced relative 
to the dust continuum maximum. They are indicated by 
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Table 1. Characteristics of the maps: telescope, line observed, half-power beamwidth, sampling step (or pixel size) in arc sec 
and in AU at the distance of the sources (150 pc), the effective resolution in pixels and in AU 



Telescope 


Line 


HPBW 


Pixel Size 


Resolution 


IRAM (30 m) 
CSO (10.4 m) 


CO (J=l-0) 
CO (J=2-l) 


22" 
28" 


7.5" (1125AU) 
16" (2 400AU) 


3 pixels (3 375 AU) 
2 pixels (4 800AU) 



Table 2. Characteristics of the mapped fields: sizes in pc, number of spectra and average signal-to-noise ratios in peak 
temperature and integrated area (The associated uncertainties are the rms scatter of the values over the maps) . 



Source 


Line 


Field 


size (pc) 


Nb. Spectra 


SNR 


SNR a 


Polaris 


12 co (J=l-0) 


0.26 


x 0.35 


3300 


11± 3 


51± 13 


Polaris 


13 CO (J=l-0) 


0.22 


x 0.30 


1650 


12± 4 


40± 12 


L1512 


12 CO (J=l-0) 


0.22 


x 0.44 


3200 


15± 3 


46± 8 


L1512 


13 CO (J=l-0) 


0.22 


x 0.36 


2560 


23± 5 


61± 13 


L1512 


12 CO (J=2-l) 


1.05 


x 1.12 


8300 


10± 4 


27± 12 




Fig. 1. Maps of the Polaris (MCLD123+24.9) field in the 
12 CO(J=1-0) (left) and 13 CO(J=1-0) (right) lines: (top) line 
integrated areas, (bottom) line centroid velocities. The offsets 
are in arcsec respectively to the (0,0) position: /// = 123.68°, 
bu — 24.93°. The two upper-left corners of the 12 CO and 
13 CO maps coincide. The dashed contour level at 4Kkms~ 1 
in the W13 map is meant to localize the region of largest 
column density and is reproduced in the Wi2qq map. The 
black contour is the 0.6 Kkms _1 -level of the CS(3-2) map 
(see Fig. 4) and is meant to localize the dense core in the rest 
of the paper. The two crosses show the positions of the HC3N 
peaks (Heithausen et al. 2002). The linear scale is shown. 



two crosses in the 13 CO map of Fig. 1. H2 densities as large 
as a few 10 5 cm~ 3 in this dense core are derived from the 
multi-line analysis of Heithausen et al. (2002) and Gcrin 
et al. (1997). Our unpublished map of CS(3-2) integrated 



"CO (.1=1 0) ] ^0O (.T=l 0) 




Fig. 2. Same as Fig. 1 for the L1512 field observed at IRAM- 
30m. Offsets are in arcsec relative to the (0,0) position: 
RA(1950)= 05h00m54.5s, Dec(1950)=32°39'00". The dashed 
contour level (3Kkms _1 in the Wi3qq map) is reproduced in 
the Wi2qq map and helps localize the regions of largest col- 
umn density. The dense core here (black contour) is identified 
with the O.lSKkms" 1 contour of the HCO + (3-2) emission (see 
Fig. 4) and is reproduced in the rest of the figures. The linear 
scale is shown. 
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Fig. 3. Maps of the 12 CO(J=2-l) line integrated area W (top) 
and line centroid velocities C (bottom) of the L1512 field ob- 
served at the CSO. In both panels, a map of the same quan- 
tities computed with the high angular resolution data of the 
IRAM-30m is inserted in the CSO map at the relevant posi- 
tion. The dense core is identified as in Fig. 4 (black contour). 
The dashed contour corresponds to Wi2qq = 2.5Kkms~ 1 
and will be used in Fig. 11. The box is the location of the 
spectra map shown in Fig. 5. The (0,0) position of the CSO 
map is RA(1950)= 05h00m54.5s, Dec(1950)=32°36'00". The 
linear scale is shown. 



emission is displayed in Fig. 4 and delineates the region 
of largest density. The contour level at 0.6 Kkms _1 of the 
CS(3-2) emission is shown in Fig. 1 to help localize the 
dense core. The 13 CO(J=1-0) contour level of 4 Kkms -1 
is also drawn on Fig. 1 to indicate the region of column 
density larger than ~4x 10 21 cm~ 2 . 

In L1512, the dense core has been mapped in CS(2- 
1) by Fuller (1989) over several arc minutes. Lee et al. 
(2001) have mapped it over a more restricted area in the 
CS(2-1) and N 2 H+ lines. A map of the HCO+(3-2) emis- 



Polaris CS (J=3-2) 




0.1 pc <5RA (") 



L1512 HCO+ (.1=3-2) 




<5RA (") 



Fig. 4. Map of integrated emission of the dense cores: IRAM- 
30m map of Polaris in the CS(3-2) line (top) and CSO map of 
L1512 in the HCO+(3-2) line (bottom). 

sion has been performed at the CSO in November 2002. 
The HCO + (3-2) line integrated emission is displayed in 
Fig. 4. As in other dense cores, the N2H+ emission is 
more concentrated than the CS emission and the N2H" 1 " 
and HCO + (3-2) emissions have very similar boundaries. 
Again, in the following, we adopt the 0.15 K km s -1 con- 
tour of the HCO + (3-2) emission as the boundary of the 
L1512 dense core. It is drawn in Figs. 2 and 3 to help 
locate the dense core. The maps of Figs. 1 to 3 confirm 
that dense cores, as traced by molecular lines such as CS, 
N2H + or HCO + at millimeter wavelengths, are embedded 
in larger structures of moderate column density, traced by 
the 13 CO lines, and further down in column density, by 
the 12 CO lines. Note that the 12 CO integrated emission 
in both fields is far from being isotropically distributed 
around the dense cores. 

The analysis presented in the next sections has been 
conducted on the 12 CO and 13 CO data sets, when avail- 
able. 12 CO is the only available tracer of molecular gas 
down to H 2 column densities as low as a few 10 20 cm~ 2 . 
Reaching these limits is mandatory in order to be able to 
make a link with statistical works performed on HI data 
(e.g. Miville-Deschenes et al. (2003); Brunt (2003)), and to 
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„',-■) 224.00 208.00 192.00 176.00 160.00 144.00 128.00 112.00 
<5Dec (") 



Fig. 5. Representative spectra of the 12 CO(J=2-l) CSO map, 
located in the box drawn in Fig. 3. On each spectrum, the 
window selected by the procedure designed to optimize the 
signal-to-noise ratio of the line integrated area is shown as 
thin vertical bars, as is the computed line centroid. The line 
temperature and velocity scales are given in the lower left cor- 
ner. 

obtain large maps avoiding star forming regions. However, 
the large opacity of the 12 CO line prevents the sampling 
of all the gas on the line of sight, most severely in regions 
close to the dense cores. The 13 CO lines are thus used 
to analyze the gas in regions of intermediate column den- 
sity between the transparent cloud edges traced by 12 CO 
only and the dense cores. The areas available for statistical 
analysis in this line are smaller than those available in the 
12 CO line because of the different abundances of the two 
isotopomers. The two lines are therefore complementary. 

4. Statistical analysis 

4.1. Computation of line centroids and their 
increments 

Centroid velocities have been computed according to 
the algorithm described in Pety (1999). The line win- 
dow used to compute the line centroid is critical (see 
Appendix A.l) and we have adjusted it locally to max- 
imize the signal-to-noise ratio of the integrated area de- 
fined as SNR a =£™Ti/v / mcr where Tj are the tempera- 
tures of the to channels within the window and a the rms 
noise level of the spectrum. As the optimal window is ex- 
pected to vary softly from one position to another, we 
smoothed the variations of the window edges with a me- 
dian, moving boxcar filter of size 5x5 pixels. The first 
moment is then computed on each spectrum using the 
optimal window found for its location. Fig. 5 illustrates 
the optimal windows found by this algorithm for a subset 
of the CSO map. This method also ensures that there is 
no significant emission out of the hnal optimal window. 
The line centroid increments are then computed over a 
lag I expressed in pixels of the map. These increments 
are computed only between spectra that (i) have a SNR a 
greater than 10.0 and (ii) have at least 4 neighbors with 



Table 3. Quantities used to normalize the <5C-PDFs for each 
field and line and for the lag I = 3 pixels: the offset applied 
to the increments to get a centered PDF (col. 3) and the stan- 
dard deviation of the SC values (col. 4). The ratio of these two 
quantities is given in col. 5. 



Field 


Line 


(5C 3 ) 
kms™ 1 




°f>c 3 
kms -1 


(sc 3 ) 

°6C 3 


Polaris 


12 C0 (J=l-0) 


8.3 xl(T 


3 


0.11 


0.07 


Polaris 


13 CO (J=l-0) 


1.5 xl(T 


2 


0.10 


0.15 


L1512 


12 CO (J=l-0) 


1.2 xl(T 


2 


0.05 


0.24 


L1512 


13 CO (J=l-0) 


1.8 xl(T 


3 


0.04 


0.04 


L1512 


12 CO (J=2-l) 


1.3 xl(T 


2 


0.09 


0.015 



SNR a > 10.0. We show in Appendix A. 2 how this thresh- 
old SNR a = 10.0 ensures that our statistical analysis of 
velocity increments is not contaminated by thermal noise. 
For a given lag I, the increments are computed for half 
the possible orientations of the lag vector because if all 
the orientations are kept for each lag, each increment in 
the map appears twice with opposite signs. In the case of 
I = 3 pixels, for instance, there are 16 different orienta- 
tions, defined by the 16 neighbors of any position lying 
within the rings of radii 2.5 and 3.5 pixels. We kept only 9 
different orientations for each position in the maps. Note 
that the smallest significant lag is I ~ 2 pixels for the 
CSO map and I — 3 pixels for the IRAM maps because 
the maps are oversampled. 

4.2. Non-Gaussian wings in the centroid increment 
PDFs. 

We have built the PDFs of the computed centroid incre- 
ments for four different lags in each field and for each line. 
Hereafter, the PDFs of centroid increments will be called 
<5C-PDFs or increment PDFs. They may also be seen as 
shear-PDFs since, as said in Section 2, the measured pro- 
jection of the velocity is orthogonal to the plane in which 
displacements are measured (i.e. the plane of the sky). 

Figs. 6 and 7 show the evolution of the (5C-PDFs with 
increasing values of the lags for each map. The PDFs are 
normalized to zero mean and unit dispersion and a nor- 
malized Gaussian PDF with the same dispersion is shown, 
for comparison, as a dotted line. The quantities used for 
this normalisation are given in Table 3 for each map and 
a lag of 3 pixels. (SC^) is the offset applied to the in- 
crements at lag I = 3 pixels to get a centered PDF and 
<Jsc 3 is the standard deviation of the increments, for the 
same lag. The ratios are given in Col. 5: they are 

all significantly smaller than unity. We therefore feel con- 
fident that the normalisation we apply does not affect the 
estimate of the <5C-PDFs computed at small lags. The er- 
ror bars drawn simply reflect the number of elements in 
each bin of the PDF. Note that the statistical significance 
of the largest increments is best for the CSO field be- 
cause the number of independent spectra available in this 
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Normalized Centroid Increments 

Fig. 6. PDFs of centroid velocity increments for increasing lag 
values I given in arcsec and pixels within each frame. The PDFs 
have been computed from a) the 12 CO (J=l-0) and b) 13 CO 
(J=l-0) spectra maps of Polaris. The centroid increment scales 
have been normalized so that (5C) = and use = 1- The 
actual velocity dispersion of each distribution cTgcu) is given in 
kms -1 in each frame. The dotted curves represent a Gaussian 
of zero mean and unit standard deviation. 

map is about nine times larger than in the IRAM maps. 
This illustrates the importance of the size of the maps 
for this kind of statistical study where departures from 
Gaussianity are searched for and occur only with proba- 
bilities close to or below. In Appendix B, we compute 
the largest statistically significant lags for each map and 
two values of the numbers of bins in the PDFs, 30 and 
8 bins. These numbers confirm that PDFs built with 30 
bins, as in Figs. 6 and 7, are statistically significant up to 
lags of ~ 6 pixels for the small fields and 12 pixels for the 
large one. They also show that the broad features in the 
PDFs, such as non-Gaussian wings, spreading over 4 bins 
or more, are significant up to lags of 12 pixels for all the 
fields. 

All the sets of increment-PDFs exhibit non-Gaussian 
wings with departure from Gaussianity being more pro- 
nounced as the lag becomes small. The effect is the most 
visible in the large scale L1512 field, because the large 
number of independent spectra in that field reduces ef- 
fects due to the poor sampling at large lags and makes 
the PDFs there the most symmetrical. 

4.3. The influence of large-scale velocity gradients 

Large scale variations of the line centroid velocities are 
visible on Figs. 1 to 3. They most likely trace large scale 
velocity gradients. The observed values are derived from 
the two extreme values of the centroid velocity C m i n and 
Cmax and the size scale over which they are measured s 
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been computed from the a) 12 CO(J=1-0) and b) 13 CO(J=1-0) 
IRAM-30m maps and c) 12 CO (J=2-l) CSO map. 

(Table 4). These numbers provide an estimate of the ob- 
served velocity gradients Vu bs ~ |C max — C m i n | /s given 
in column 5 of Table 4. We note here that the values de- 
rived from two different lines in the same field are not 
necessarily the same because the lines do not quite sam- 
ple the same gas, as discussed above. Unlike in other stud- 
ies (Grosdidier et al. 2001; Miesch et al. 1999), we did not 
remove the contribution of large scale velocity gradients 
in our centroid increment computations. In this section, 
we explain why. 

In the analysis of turbulence in HII regions, the large 
scale velocity gradients are related to the expansion of 
the warm ionized gas, driven by the pressure gradient of 
the HII region relative to the surrounding medium. It is 
thus justified to remove large scale gradients, prior to any 
statistical analysis of the turbulence within HII regions, as 
done by Grosdidier et al. (2001), because the HII region 
expansion is an ordered large scale motion and is not part 
of the inertial range of the turbulent cascade. 

Miesch et al. (1999), who analyze the turbulence in 
molecular clouds, have different goals from ours. Their 
goal is to extract an homogeneous and isotropic subset 
of interstellar turbulence and compare it to laboratory 
experiments or direct numerical simulations. Instead, we 
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Table 4. Influence of large scale velocity gradients. Difference between the observed extrema of line centroids in each map 
(col. 3), the associated scale s (col. 4) and the corresponding observed large scale velocity gradients (col. 5), the contribution 
of this large scale gradient to a lag I — 3 pixels (col. 6), the ratio of this contribution to the internal velocity dispersion at the 
same scale (col. 7), the exponent of the scaling of the standard deviation of the centroid increments with the lag at which they 
are computed (col. 8), the velocity gradients at scale s inferred from this scaling and ascribed to turbulence (col. 9). 
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Fig. 8. Scaling of the standard deviations of centroid incre- 
ments with the lag I used to compute the increments for the five 
fields and lines. Large symbols are used for the lags which give 
statistically meaningful results as explained in Appendix B. 
The exponents £2 listed in Table 4 have been determined for 
the range of lags which is statistically meaningful. 

search for regions of enhanced shear at small scales to lo- 
calize bursts of dissipation. We show below that we are 
not limited in our analysis by large scale anisotropy. We 
first argue, following the results of numerical simulations 
of Burkert & Bodcnheimcr (2000) that turbulence itself 
with its steep power spectrum (i.e. most of the power is 
in the large scales) generates velocity gradients which may 
be interpreted as rotation at any scale. We show in Fig. 8 
and Table 4 that this is indeed the case. The scaling of the 
standard deviation of the line centroid increments asci 
with lag / is shown in Fig. 8. Up to lags of ~ 100 arc sec 
or about 12 pixels, asci a The values of (2 are given 
in Table 4. This scaling provides that of the turbulent ve- 
locity shear V«t U rb l^ 2 ^ 1 with I. The turbulent velocity 
shear, Vw tU rb, inferred from this scaling law at the scale 
s of the largest centroid differences is given in col. 9 of 
Table 4. Comparison of columns 5 and 9 shows that in 
all cases: Vv bs & Vv tU rb- The large scale observed veloc- 
ity gradients are therefore smaller than or comparable to 



those expected from the turbulent shear at the same scale, 
suggesting that they are part of the turbulent dynamics. 

Last, we show that the large-scale gradients do not af- 
fect the small scale statistics which are of importance to 
the present analysis. The contribution of the large-scale 
gradients to the centroid velocity increments computed as 
Sv 3 — Vuobs x I for I — 3 pixels are given in Table 4 for 
each map (col. 6). These values are then compared to the 
standard deviation of the 5C-PDFs distribution (col. 7) . 
In all cases, the ratio 8v^/asc 3 is smaller than or equal 
to unity. The contribution of the large-scale gradients to 
the centroid increments at small scale is therefore much 
smaller than the increments populating the non-Gaussian 
wings which extend up to 4 or 5 age (see Figs. 6 and 7). 
Furthermore, the large scale gradient has a well-defined di- 
rection and therefore preferentially affects the increments 
computed along the same direction i.e. one increment out 
of 9 in the case of I — 3 pixels (out of more for larger 
lags). The contribution of the large scale velocity gradient 
to increments computed over lags of different orientations 
is further reduced by the appropriate cosine. For these 
reasons, we are confident that keeping the large scale ve- 
locity gradients in the maps does not significantly distort 
the (5C-PDFs of Figs. 6 and 7. 

In summary, we did not remove large-scale velocity 
gradients prior to our statistical analysis for two reasons: 
(1) there is evidence for the large-scale gradients to be 
part of the turbulent dynamics we analyze and (2) their 
value is such that their contamination of the small scale 
statistics on which we focus here, is negligible. 

5. Spatial distribution of the non-Gaussian 

centroid increments: a new kind of small scale 
structure 

5.1. The positions populating the non-Gaussian wings 
of the increments PDFs are not randomly 
distributed: they form elongated structures 

We have seen in the previous section that, in most cases, 
the departure of the increment PDFs from a Gaussian dis- 
tribution is as large as the lag is small. We are here inter- 
ested in tracing the locations of the positions for which the 
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Table 5. Characteristics of the maps of averaged centroid increments for each field and line: thresholds used to define the 
subsets of positions populating the non-Gaussian wings of <5C-PDFs (col. 3), peak increment (col. 4) and background (col. 5) 
values, square of the contrast of the peak to background value (col. 6), average of the line profile velocity dispersion over the 
field (col. 7), surface filling factor of regions where |<5C| > \8Co\ (col. 8) and fraction of dissipation in these regions (col. 9, see 
Section 7). The last column is the average linewidth in each field (see Section 6.3). 
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Fig. 9. Top panels: Maps of the averaged centroid increments 
|<5C| in Polaris computed with a lag of 3 pixels from the 
12 CO(J=1-0) lines (left) and the 13 CO(J=1-0) lines (right). 
The solid contour localizes the dense core (see Fig. 1). The 
bottom panels are sketches meant to emphasize the relative 
orientation of the region of large column density traced by 
13 CO (grey area), the dense core (dashed contour) and regions 
of large centroid increments where |<SC| > \6Co\ (black con- 
tours). The thresholds \SCo\ are those listed in Table 5. The 
letters are the labels used in the text. 

small scale increment values build up the non-Gaussian 
wings. We display in Figs. 9 through 11 the spatial distri- 
bution of the increments computed for lags of three pixels. 
The lag value I = 3 pixels is the smallest significant lag in 
the IRAM maps because the maps are Nyquist sampled or 
better: adjacent pixels are therefore not independent (see 
Table 1). 

The spatial distributions of the centroid increments 
are shown as azimuthal averages over the 16 possible di- 
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Fig. 10. Same as Fig. 9 for the L1512 field observed at IRAM. 



rections because it allows a more compact presentation 
of the results. In the rest of the paper, we note \SC\ the 
azimuthal average of the modulus of all the oriented cen- 
troid increments computed at a given position, for lags of 
three pixels. The thresholds, |<5C 1 , given in Table 5, has 
been chosen to isolate the regions where more than half of 
the oriented increments belong to the non-Gaussian wings 
of the (5C-PDFs. These thresholds are close to the stan- 
dard deviations <rsc 3 of the increment-PDFs for I = 3 
pixels. Note that the averaging produces increment values 
smaller by factors up to a few than the original values used 
to build the PDFs. This effect is illustrated in the cuts of 
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Fig. 11. Same as Fig. 9 for the large scale L1512 field. The 
dot-dashed contour is the same as in Fig. 3. 



Figs. 13-16 where the variations of the centroid and of the 
centroid increments before and after azimuthal averaging 
are shown as a function of position. 

In Figs. 9-11, the uniform black areas correspond to 
positions where increments are not computed because at 
least one of two spectra used either has a SNR a value 
below 10.0 or is too isolated (i.e. with less that 4 con- 
tiguous neighbors with SNR a > 10.0). The black contours 
drawn after the CS(3-2) or HCO + (3-2) maps help locate 
the dense cores. 

In the sketches of Figs. 9-11, the black contours la- 
belled by letters follow the thresholds \6Cq\. These con- 
tours are therefore intended to trace the structures delin- 
eated by the positions where the original, non averaged, 
increments are larger than a few asc 3 an d belong to the 
non-Gaussian wings of the <5C-PDFs. 

Two characteristics of these structures may be derived 
from Figs. 9-11. Firstly, the spatial distribution of the 
largest averaged increments is not random. They are con- 
centrated in specific regions and appear to form elongated 




<5RA (" 



Fig. 12. Maps of average centroid increments of the same 
region of the L1512 field computed from the IRAM-30m 
12 CO(J=1-0) map (left) and the CSO 12 CO(J=2-l) map 
(right). Note the remarkable quantitative agreement between 
the two maps. 



structures, almost straight in several cases, and covering 
only a small fraction of the maps. Secondly, the \6C\ values 
in these structures are well above the background fluctua- 
tions of the increments in the maps. Table 5 gives the back- 
ground values |<5Cb g | computed as the median of all the 
values smaller than the thresholds \5Co\ and the largest 
values of increments |5C max | in the maps. It shows that 
the largest values are between 4 and 9 times above the 
background values, far above statistical fluctuations. 

Last, it is remarkable that similar spatial distributions 
are found for the large increments computed in the same 
field with two independent data sets i.e. from two differ- 
ent telescopes and two different lines, 12 CO (J=l-0) and 
(J=2-l). Fig. 12 allows a comparison of the patterns found 
for the centroid increments in the L1512 field from the 
IRAM-30m and CSO maps. The morphological and quan- 
titative agreement between the two maps fully supports 
the reality of these structures. In Section 6, we discuss 
why they do trace the shear of the velocity field. 

The smallest lags of three pixels introduce an effective 
resolution, ascribing a size of three pixels to unresolved 
variations of the line centroid. But some of the structures 
traced by the largest increments are resolved by the ob- 
servations. This last characteristic is also seen in the cuts 
of Figs. 13-16. In the RA cut across the L1512 field, the 
peak of increment visible at RA offset <~ 200" is due to 
the steepening of the variation of the line centroid over 
6 consecutive pixels, from offset <~ 300" to 200" or more 
than three CSO bcamwidths. The half-power width of this 
feature (above the background value), ~ 0.05 pc, is there- 
fore resolved by the procedure of computing increments 
over lags of 3 pixels. We illustrate with the other cuts 
that several of these structures are actually resolved, with 
thicknesses of the order of 0.05 pc, up to 0.08 pc. 
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Fig. 13. Cuts in the 12 CO line of the Polaris field in RA, aver- Fig. 15. Same as Fig. 13 for the large scale L1512 field. The 

aged over the indicated Dec range. The histograms are the cuts letters correspond to those of Fig. 11. 

of increment centroids: raw values in the top panel, averaged 

in azimuth as explained in Section 5.1 in the bottom panel. 

The other cuts are the centroid velocities (top) and integrated 

area (bottom). The letters correspond to those of Fig. 9. 
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Fig. 14. Same as Fig. 13 for the 13 CO line of the Polaris field. 
Here the cuts are in RA, averaged over the indicated Dec range. 



Fig. 16. Same as Fig. 14 for the large scale L1512 field. 
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5.2. Comparison with the results of Miesch et al. 
(1999) 

The non-random distribution of the largest centroid veloc- 
ity increments seems at odd with the spotty distribution 
found by Miesch et al. (1999) in their maps of centroid ve- 
locity differences of several star forming regions. However, 
we consider that our results agree with theirs for the fol- 
lowing reason. Their method, as explained in Section 4.3, 
is different from ours because they first remove a smoothed 
value from the line centroid values in their maps. The cen- 
troid velocity fluctuations they obtain with this procedure 
are already of the same nature as a velocity increment, as 
shown by Lis et al. (1996). It is therefore meaningful to 
compare their maps of centroid fluctuations with our maps 
of centroid increments. Their maps of centroid velocity 
fluctuations (their Figs. 4 and 5) do indeed exhibit spa- 
tial structures. The largest values appear in several cases 
to form narrow and elongated structures similar to those 
we find. The spotty distribution is found only once they 
compute the increments of these former maps of velocity 
fluctuations. What they call centroid velocity differences 
cannot be compared to what we compute. 

A more detailed comparison of the two sets of results 
would not be meaningful for several reasons. The signal- 
to-noise ratios of our data given in Table 2 are signifi- 
cantly better than those of Miesch et al. (1999). Their 
SNR values range between 4.5 and 8.9, except for HH83, 
which has SNR=19.7 but is their smallest field. The na- 
ture of the fields is also different. While we analyze only 
quiescent fields in the present paper, Miesch et al. have 
analyzed active star forming regions where many shocks, 
driven by the young stars, interact and probably generate 
considerable small scale structure in the patterns of the 
vorticity. Last, most of the fields being at distances larger 
that the nearby fields studied here, the angular character- 
istic scales of the structures (if any) are expected to be 
smaller. 

5.3. The regions of largest centroid increments are not 
coinciding with density nor column density peaks 

The sketches of Figs. 9 to 11 allow a comparison of the 
structures of largest increments with the concentrations 
of matter, whether they are large column density or large 
density structures or both. The dotted contours delin- 
eate the dense cores and the grey areas those where the 
13 CO(J=1-0) integrated emission is larger than 3 and 4 
Kkms -1 , for L1512 and Polaris respectively, which corre- 
sponds to N(R 2 ) > a few 10 21 cm~ 2 , about 10 times the 
average column density of the large scale environment. 
The main result visible in those figures is that none of the 
regions of largest increments coincides, even in projection, 
with either a dense core (i.e. a density peak) or a column 
density peak. 

This result is better seen in the cuts of Figs. 13-16: 
the large variations at small scale of the line centroid, at 
the origin of the large increment values, do not coincide 
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Fig. 17. Scatter plots of the integrated emission versus the 
centroid increments for three data sets analyzed. 



with column density extrema. Such a coincidence would 
be expected if the line centroid variations were due to the 
overlap on the line of sight of unrelated gas components. 
In that case, an extrcmum of column density would coin- 
cide with the region of large centroid increment, which is 
that over which the two components arc intercepted, pro- 
ducing a modified (locally broader) line profile. This lack 
of correlation is also seen in the scatter plots of Fig. 17 
where the largest increments are in no case associated with 
extrema of column density. Conversely, the largest incre- 
ments tend to occur at the positions were the integrated 
areas are the smallest. 

5.4. Link between the orientation of these structures 
and that of density structures 

The previous section shows that the structures traced by 
the regions of largest centroid increments at small scale 
do not coincide (in projection) with the dense cores nor 
the column density extrema in the maps. 
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Nonetheless, this statement does not mean that these 
small scale structures are unrelated to the dense cores. 
The orientation of the most elongated structures tend to 
be related with that of the dense cores. In the Polaris field, 
the structure labelled C in the 12 CO map, appears as a 
prolongation of the bright 13 CO structure, which seems 
to be an extension of the dense core itself (see sketches 
of Fig. 9). The structures B and D' are adjacent and par- 
allel to the edge of the dense core, as is seen also in the 
cuts of Fig. 14. In the L1512 field (Fig. 11), similar pat- 
terns are visible. The two southern structures A and A' 
frame in projection the region of bright 12 CO lines (dot- 
dashed contour), and the structure B is elongated in the 
same direction as that of the filament of matter seen in 
13 CO and 12 CO and the orientation of the dense core it- 
self as mapped in HCO + . Since the velocity is continuous 
between the structures of large increments and those of 
matter (see maps of line centroids in Figs. 1-3), it is most 
likely that they are connected in the 3-dimcnsional space. 

5.5. Volume filling factor of the largest increment 
structures 

For each field, in addition to the maximum value of the 
increments in the map, |5C max | and the background value 
|5Cbg| defined in Section 5.1, Table 5 also gathers the frac- 
tion f 8 of pixels in the map where the centroid increments 
are larger than the thresholds |5Cb|. f s may be under- 
stood as their surface filling factor in the maps. From its 
values, one derives volume filling factors of the regions of 
largest increments of at most a few %. Assuming that the 
depth of the structures along the line of sight is compa- 
rable to their projected thickness, d, and that the extent 
of the cloud along the line of sight, L, is provided by the 
large scale maps (Falgarone et al. 1998), f v ~ f s d/L. For 
d <~ 0.05 pc and L <~ 1 pc, /„ <~ a few % or less. 

6. Discussion 

The subset of spectra shown in Fig. 5 encompasses one of 
the regions where the centroid increments are the largest. 
Fig. 5 shows that the optimal window found to compute 
the line centroids varies only by very small amounts from 
one spectrum to the next. This illustrates that the mea- 
sured centroid increments are due to small variations in 
the shapes of the line profiles from one line of sight to the 
next. These variations may have several origins, unrelated 
to the structure of the velocity field: optical depth effects, 
gas temperature or density fluctuations. We show in the 
following section that our interpretation of these varia- 
tions as tracers of small scale velocity shears is justified. 

6.1. Optical depth effects 

The optical depth of the lines used introduces a possible 
bias in the sampling of the velocity field performed by 
the spectra. Indeed, optically thick lines tend to sample 
preferentially foreground material. 



However, the maps of centroid increments built from 
the 13 CO and 12 CO lines, although not exactly the same, 
are reminiscent of each other, suggesting that the two 
lines sample about the same gas (see Figs 9-10). This is 
also suggested by the similarity of the maps obtained in 
L1512 with two different 12 CO transitions of different op- 
tical depths (Fig. 12). These results support the concept 
of low effective optical depth of the 12 CO lines discussed 
in Martin et al. (1984). Falgarone et al. (1998) also in- 
ferred a low effective optical depth in the 12 CO lines of 
the fields analyzed here (except in the projected area of 
the dense cores), on the basis of the uniformity of the 
excitation temperature. Note that a low effective optical 
depth of the 12 CO lines is to be expected in these fields (as 
opposed to brighter molecular clouds) because a large ve- 
locity dispersion of the gas, leading to a large escape prob- 
ability of the photons in velocity space, is combined with 
low gas column densities. The same similarity is found 
between maps of centroid increments built from the 13 CO 
and C 18 spectra of the environment of the L1689B dense 
core (Pety et al. in preparation), suggesting that optical 
depth effects are not seriously affecting the centroid incre- 
ment determinations. 

This result illustrates a major strength of the line cen- 
troid analysis. On the one hand, centroids are velocity mo- 
ments and as such provide more weight to velocities far 
from the line centroids, i.e. the line wings. The wings are 
the parts of the line profiles the least affected by optical 
depth effects. On the other hand, our method which limits 
the range of velocities used to compute the centroid (see 
Section 4.1), limits the noise contribution, particularly in 
the line wings. We thus argue that there is a velocity range 
in spectrally well-sampled line profiles, neither too far in 
the line wings (to avoid noise artefacts) nor too close to 
the line centroid (to avoid optical depth artefacts), which 
carries enough information on the velocity field to be used 
on statistical grounds. 

Last, self-absorption is also unlikely in those low 
brightness core environments because the excitation tem- 
perature of the 12 CO and 13 CO lines is remarkably uni- 
form as was shown by the scatter plots of the (J=l-0) 
to (J=2-l) 12 CO and 13 CO lines (Falgarone et al. 1998). 
For self-absorption to occur the optical depth must ex- 
ceed unity and the foreground material must have a lower 
excitation temperature. 

6.2. Role of temperature fluctuations 

It might be argued that the variations of line centroids 
traced by the 12 CO lines are due only to variations of the 
gas kinetic temperature. This might be the case in bright 
sources where the 12 CO lines are thermalized. But it has 
to be appreciated that in the fields under study here the 
12 CO lines are weak (less than 5 K) while the gas is poorly 
shielded from the ambient UV field (H2 column densities 
less than a few 10 21 cm~ 2 ) and thus likely to be much 
warmer than 5K. It means that either the lines are sub- 
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thermally excited, in which case their excitation is mostly 
radiative, or that they are thermally excited in small struc- 
tures with a large beam dilution factor. In either case, the 
line profile variations are not determined by variations of 
the gas kinetic temperature but by those of the radiative 
excitation and/or beam filling factor, which both depend 
on the velocity field (Falgarone et al. 1998). 

6.3. Role of density fluctuations unrelated to the 
velocity field 

In this section, we discuss the influence of column density 
or density variations, independent of the velocity field, on 
the velocity increments that we measure. Under the op- 
tically thin hypothesis, dN(u) oc T(u) du where dN(u) 
is the column density of gas at the projected velocity u 
within du. The centroid velocity defined in Eq. (1) is there- 
fore also 



C = 



J udN(u) 



N 



with N 



= J dN(u). 



The centroid increment between lines of sight 1 and 2 can 
thus be written 



J uidJVi(ui) 



/ u 2 dN 2 (u 2 ) 

w 2 



(4) 



From this expression, it is obvious that 8C\- 2 = 
if the two lines of sight differ only by their total col- 
umn density, the velocity structure being the same i.e. 
N 2 (u 2 ) oc N\(u\) for all velocities. In this case the two 
integrated line profiles are simply homothetic. 

Let us consider another simple case in which the total 
column density N along the two lines of sight is approx- 
imately the same, the only difference being that a frac- 
tion of the gas at projected velocity u\ on the first line 
of sight is at u 2 on the other. This is possible, for in- 
stance, in the perspective of matter distributed in a large 
number of very small structures, possibly virialized in the 
potential well of the cloud. Then 5G\- 2 = (u\ —u 2 ) SN/N 
where SN is the column density of the gas which is at 
different projected velocities on each line of sight. We es- 
timate below the value of SN/N required to produce cen- 
troid increments of the order of |<5C max | in each field. We 
assume that the velocity difference u\ — u 2 is of the or- 
der of magnitude of the average linewidth in each field 
(Aw) given in Table 5, Ui — u 2 ~ (Aw). Therefore an 
estimate of the required column density fluctuations is: 
SN/N <~ |<5C max | / (Aw). The latter ratio can be computed 
from the values in Table 5 and is equal to 1.25, 0.5 and 
1.05 for the three fields of Fig. 18, respectively. Now, if we 
assume that the relative variations of the line integrated 
areas trace those of gas column density, we can check how 
the observed fluctuations of column density compare with 
the above values. This may be seen in the scatter plots of 
the averaged increments versus the relative variations of 
the line integrated area, |<L4| /A (Fig. 18). These scatter 
plots are shown only for the three maps where the statis- 
tics are the least affected by the dense core. The values of 
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Fig. 18. Scatter plots of the azimuthally averaged increments, 
\SC\, versus the relative variations of the line integrated area, 
|<L4|/A, for three data sets. The horizontal lines show the 
| S Co | thresholds separating the Gaussian kernel from the non- 
Gaussian wings in the <5C-PDFs. The vertical lines are an es- 
timation of the required column density relative variations to 
produce centroid increments of the order of |<5C m ax| in each 
field. 



|<5Cmax| /(Aw) and of the threshold |<JC | are indicated for 
each field. The scatter plots show that, except for a few 
lines of sight of the Polaris 13 CO map, the observed rel- 
ative variations of column density (or |<L4| /A) at the po- 
sitions of increments populating the non-Gaussian wings 
(i.e. increments above the thresholds \8Cq\) are smaller 
than the minimum value estimated above. In the large 
scale L1512 field and the Polaris field as seen in 12 CO, 
it is therefore unlikely that the observed centroid incre- 
ments trace small scale column density fluctuations due 
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for instance to randomly distributed clumps. This is also 
supported by the fact that the structures traced by the 
large centroid increments form spatially coherent patterns 
and are not randomly distributed as would be expected for 
randomly moving clumps. We thus conclude that the con- 
tribution of small scale random column density variations 
to the observed centroid increments is not dominant. 

These small scale column density fluctuations could 
trace density fluctuations. They occur over lags of 3 pix- 
els or I ~ 5000 AU (see Table 1). The corresponding den- 
sity enhancements are therefore expected to occur over the 
same small scale I. Thus, Sn/n <~ (\5A\ /A) (L/l) where L 
is the depth of the cloud. For \SA\ /A = 0.5 and L = 1 
pc, 5n/n = 25. These might be supersonic shocks with 
Mach numbers of about 5, in which case the density and 
velocity fields are closely correlated. However, such large 
density enhancements would produce variations of the 
CO(J=2-1)/CO(J=1-0) line ratio in the weak regions of 
the fields, variations which are not observed (Falgarone 
et al. 1998). 

The tenuous variations of the line profiles at the origin 
of the large centroid increments at small scale are therefore 
most likely ascribed to large velocity shears at small scale. 

7. Possible link with the dissipation of turbulence 

Our statistical analysis of the environment of two low 
mass dense cores leads to two main results: (i) the PDFs 
of line centroid increments (or velocity shears) have non- 
Gaussian wings as more pronounced as the lag over which 
they are computed is small, and (ii) the spatial distribu- 
tion of the largest shears is not random and reveals often 
resolved and elongated structures not coinciding with den- 
sity nor column density peaks. These are predominantly 
velocity structures. 

The first property is a characteristic shared by incom- 
pressible turbulence, compressible turbulence and com- 
pressible, magnetized turbulence and is a signature of the 
intcrmittency of the velocity or current density fields. Our 
results therefore suggest that the observed distributions of 
large shears may trace regions of intense vorticity reminis- 
cent of those responsible for the intermittency of turbu- 
lence. Note that these regions of intense vorticity are cm- 
bedded in the bulk of the flow and that the large centroid 
increments are only tracers of their presence in the flow. 
In other words, a small scale structure of large vorticity, 
embedded in the bulk of the flow, can be at the origin of a 
small variation of the line centroids, along the boundary of 
two large scale regions of more uniform centroid velocity, 
i.e. the large eddies. 

We have argued that the large shears are unlikely to 
trace shocks but we cannot rule out the fact that the vor- 
ticity extrema are fossil structures of shocks, already dis- 
rupted. 

If the largest shears that we have detected trace local 
enhancements of the vorticity and therefore dissipation, 
as discussed in the Introduction, an estimate of the lo- 
cal enhancement of the dissipation rate in those regions 



is provided by the ratio (<5C max /<5Cb g ) 2 (Table 5) be- 
cause the dissipation rate follows the square of the vor- 
ticity (or of the velocity shear) (Landau & Lifshitz, 1987). 
The local enhancements are large, between 15 and 80 
(Table 5). It is interesting to note that the smallest values 
of (<5C max /i5Cbg) 2 are obtained in the close vicinity of the 
L1512 dense core (IRAM field) and the largest in Polaris 
and in the large scale environment of L1512 (CSO field). 
This may reflect the fact that turbulent dissipation is a 
less violent process within a dense core which is already 
formed. We also note that even the weakest structures in 
the maps of \SC\ have a contrast larger than 2 above the 
background and therefore correspond to a dissipation rate 
locally 4 times larger than in the bulk of the field (see 
for instance the weak structures, labelled B and C, in the 
North of the IRAM L1512 12 CO map). 

Last, the quantity f e = T,\ S c\>\5C \ \ SC \ VEaii \ 6C \ 2 
(Table 5) may be understood as the fraction of the dissipa- 
tion taking place in the regions of large increments (those 
populating the non-Gaussian wings). The underlying as- 
sumption is that the ratio of the actual dissipation rate to 
\SC\ 2 is the same in the regions where the increments 
populate the Gaussian cores of the PDFs (the bulk of 
the fields) as in those populating the non-Gaussian wings. 
Under this assumption, the values of Table 5 show that a 
significant fraction (more than half) of the total dissipa- 
tion takes place in the ensemble of small scale structures 
filling a volume close to a few % of the cloud volume, as 
estimated above. 

8. Conclusion 

We have found a new kind of small scale structure in the 
environment of low mass starless dense cores. 

They emerge in the maps of line centroid increments 
between adjacent spectra as regions between 4 and 9 times 
brighter than the background fluctuations. They are es- 
sentially velocity structures because they are not coincid- 
ing with any detected sufficient increase in the column 
density or density. They are small scale, often elongated, 
structures. A few are resolved by the observations with 
thicknesses up to 0.08 pc. Most of them are not resolved 
i.e. smaller than 0.02 pc, the effective resolution due to 
the smallest lag available to compute the increments. 

The measured projected velocity being orthogonal to 
the plane of the sky, in which lags are measured, the cen- 
troid velocity increments trace one component of the shear 
of the velocity field. The shear-PDFs exhibit non-Gaussian 
wings which are the most prominent for the smallest lag. 
For this reason, we ascribe these structures to turbulence 
by analogy with the behavior of such PDFs in labora- 
tory flows or numerical simulations of compressible tur- 
bulence. The new structures delineate the locus of posi- 
tions where the increments are the largest and populate 
the non-Gaussian wings. 

These structures therefore possibly trace regions of en- 
hanced dissipation of turbulence in the observed fields. 
The nature of the dissipation process is not determined 
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by the observations: it could be cither viscous dissipation, 
ion-neutral collisions or enhanced current densities. These 
structures are unlikely to trace shocks, but they may trace 
the remanent vorticity generated by shocks already dis- 
rupted. In this picture, a significant fraction of the tur- 
bulent energy present in the field would be dissipating in 
those structures filling less than a few % of the cloud vol- 
ume. 

Last, the structures display some connexion with the 
tracers of dense gas, such as shared orientations and ve- 
locity and space pattern continuity. Before we can make a 
link between such structures and the formation of dense 
cores, we need to investigate control fields i.e. molecular 
clouds without dense cores (Hily-Blant et al. in prepara- 
tion) and with moderate star formation activity (Pety et 
al., in preparation). 
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Appendix A: Role of the thermal noise in the 
statistical analysis 

Thermal noise in the spectra affects the estimation of the line 
ccntroid velocity. It thus could be that thermal noise largely 
contributes to the non-Gaussian wings of the 5C-PDFs. This 
influence of thermal noise is relatively more important at small 
lags where the increments are the smallest because of the scal- 
ing laws of turbulence. In section A.l, we show that most of 
the large increment values are significantly above the thermal 
noise contribution. And in section A. 2, we present results of nu- 
merical simulations performed to answer two questions: (i) can 
thermal noise create non-Gaussian wings on intrinsically (orig- 
inally) Gaussian <5C-PDFs and (ii) how does thermal noise af- 
fect the shape of an <5C-PDF which has intrinsic non-Gaussian 
wings? 



A.l. Increment SNR 

The uncertainty on the centroid velocity of a spectrum scales 
as (Pety 1999): 



AC oc VF/SNRa 



(A.l) 



where W is the width of the optimal window (see Section 4.1). 
The uncertainty on the centroid velocity increments, A(5C), 
can be approximated as the quadratic sum of the uncertainties 
on the centroid of each spectrum in the pair. Fig. A.l shows the 
scatter plots of the centroid velocity increments signal-to-noise 
ratios (SNR^c, defined as the ratio of SC to its uncertainty) 
versus the centroid velocity increments for the five maps. Here 
all the directions of the pairs are considered. 

The butterfly shape of these scatter plots is due to the 
fact that the centroid uncertainties are bounded, essentially by 
the two extreme values of SNR a , since W does not vary much 
in a map. We impose the lowest SNR a value to be 10.0 (see 
Appendix A. 2). The upper value reflects the finite noise level 
and the peak line intensity in each map. The <5C-uncertainties, 
A(SC), defined as above, are therefore also bounded. We call 
these boundaries A and B so that 1/A < A(SC) < 1/B. We 
then obtain 



B * \SC\ < SNRsc 



SC 



A(SC) 



< A* \SC\. 



(A.2) 



The A and B boundaries thus appear as the slopes of the lines 
that limit the scatter plots of Fig. A.l. 

These scatter plots also show that the structures delineated 
by the largest centroid increments and shown in Figs. 9 to 11 
have large SNRs. Clearly these structures are far above the 
noise level, as was already suggested by the good agreement 
between the IRAM and CSO results for the L1512 core (see 
Fig. 12 and section 5.1). 

A.2. Simulations 

We computed two spectra maps made of Gaussian lines with 
the same line areas and widths at each position as in the L1512 
large scale field. The first one also has the same spatial distri- 
bution of centroid velocities as the L1512 field. For the second 
one, the spatial distribution of centroid velocities has been 
computed from a random, Gaussian data cube whose power 
spectrum is a power law of -5/3 exponent. It has been checked 
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Fig. A.l. Scatter plots of the SC signal-to-noise ratio (defined 
as the ratio of SC to its uncertainty) versus the 5C-values for 
the five data sets analyzed. These scatter plots show that the 
structures delineated by the largest centroid increments and 
shown in Figs. 9 to 11 have large SNRs. 
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Normalized Centroid Increments 
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Fig. A. 2. Results of numerical simulations showing the evo- 
lution of the shape of the normalized <5C-PDFs as thermal 
noise is progressively added to maps of noiseless Gaussian line 
profiles. The two maps have the same spatial distribution of 
line areas and widths as the real data of the L1512 large scale 
field. Their difference lies in the distribution of line centroids: 
it is Gaussian (top panels) or identical to the real data (bot- 
tom panels) (see text). Two different SNR a thresholds (i.e. 3 
and 10) are used in the computation of the centroid velocity 
increments. For the Gaussian velocity field, thermal noise can 
create non-Gaussian wings to the (5C-PDFs, only when the 
SNR a threshold is 3 and for the largest noise. In comparison, 
the evolution in the bottom panels shows that the L1512 veloc- 
ity field has intrinsic non-Gaussian statistical properties. The 
Gaussian PDF in the case of largest noise and SNR a is due to 
poor sampling (see text). 



that this cube does not show intermittent behavior, i.e. SC- 
PDF are Gaussian at all scales (Pety 1999). These two spectra 
maps are respectively referred to as the L1512 velocity field 
and the Gaussian velocity field. 

To mimic thermal noise, we added random numbers sam- 
pled from a Gaussian law of zero mean and a given standard 
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Fig. A. 3. Evolution of the number of spectra whose SNR a 
value is above a given threshold (3 for the dotted line and 
10 for the plain line) as a function of the noise level. 



deviation, to the intensities of all channels. We checked that 
the average noise of the spectra in the map is the value of this 
standard deviation. The original L1512 CSO field has quite 
a uniform noise distribution of average value 0.2 K. We thus 
fixed the standard deviation of Gaussian noise to values be- 
tween 0.02 and 0.8 K to explore a fair range of SNR a (the 
spectra area being kept fixed). We then used our algorithm 
to compute the centroid map from the spectra map (see sec- 
tion 4.1). We used two different SNR a thresholds (i.e. 3 and 10) 
to thoroughly explore thermal noise effects. All the <5C-PDFs 
have been computed for a lag of 3 pixels, and normalized to 
zero mean and unit dispersion. 

Fig. A. 2 shows normalized <5C-PDFs for 3 different val- 
ues of noise levels: 0.002 K equivalent to the noise-free case, 
0.2 K equal to the value of the real L1512 data and 0.8 K, the 
largest value. For the Gaussian velocity field, thermal noise can 
create non-Gaussian wings to the <5C-PDFs. However, those 
non-Gaussian wings appear only when the SNR a threshold is 
3 and for the largest noise values. This implies that thermal 
noise cannot be responsible for the non-Gaussian wings seen 
in the 5C-PDFs computed on the real L1512 data: the L1512 
velocity field has intrinsic non-Gaussian statistical properties. 

For the L1512 Velocity Field (see Fig. A. 2), thermal noise 
does not change the normalized 5C-PDF shape for typical 
noise values (i.e. 0.2 K). <5C-PDF shapes change only for the 
highest noise values. In particular, the (5C-PDF can become 
Gaussian when the SNR a threshold is 10. This is just a sam- 
pling effect. Indeed, Fig. A. 3 shows the number of points in the 
spectra map whose SNR a values are above the threshold, as a 
function of the map noise level. It can be seen that the number 
of spectra with SNR a > 10.0 (solid line) decreases much faster 
than the number of spectra with SNR a > 3.0 (dotted line). 
Due to the anticorrelation between high area and high |<5C| 
regions (see Fig. 17.a), the remaining points naturally have a 
velocity field with a Gaussian statistic. 



Appendix B: Sampling effects 

Statistical estimators are significant only if they are measured 
over a large enough number of realizations. The number of re- 
alizations of the increments at a given lag, increases as the lag 
decreases. We compute below the largest statistically signifi- 
cant lag in each available map. We request that the largest 
lag provides at least a minimum number of 10 independent 
elements in each bin of the 5C-PDFs. 
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Table B.l. Largest significant lags for the <5C-PDFs in each 
field. The field sizes are N and M in pixels and the largest 
significant lags are given (in pixels) for JC-PDFs built with 8 
(col. 6) and 30 (col. 7) bins. 



Source 


Line 


N 


M 


^max (8) 


/max (30) 


Polaris 


rz CO (J=l-0) 


48 


64 


12 


6 


Polaris 


13 CO (J=l-0) 


40 


56 


12 


6 


L1512 


12 CO (J=l-0) 


40 


80 


13 


7 


L1512 


13 CO (J=l-0) 


56 


40 


12 


6 


L1512 


12 CO (J=2-l) 


90 


96 


21 


11 



We call Af s ig n the minimum number of independent mea- 
surements required to build a statistically meaningful PDF. 
From the above condition, iV s i gn is at least 10 times the num- 
ber of bins iVbin in the PDF. The PDFs shown in Figs. 6 and 7 
are built with 30 bins therefore N sign ^ 300. 

To compute -/V s i gn , we consider that two increments are 
independent if i) the two pairs of points have at most one 
point in common and ii) the length between the two points 
which differ is at least equal to the lag I used to compute the 
increments. This makes the reasonable assumption that in the 
inertial range of turbulence, variations of the velocity at scale 
I decorrelate over the distance I. We therefore divide the map 
in cells of size I 2 . In each cell, we find four independent pairs 
of points: the two diagonals and half the four cell edges (since 
each edge belongs to two contiguous cells). Hence, there are 
about 4NM/1 2 independent increments in a map of N x M 
pixels. Thus, 

iV sign ~ 4^ (B.l) 

from which we derive the largest lag, i max , for which the PDF 
is statistically significant, as a function of the number of bins 
in the PDF: 



^ (jv -W^- (R2) 

Table B.l gives numerical values of the largest reliable lag for 
the two cases of <5C-PDFs built with 30 and 8 bins. We see, as 
expected, that the most statistically significant PDFs are those 
obtained for the large scale field of L1512, and that at least for 
the two smallest lags (i.e. 3 and 6 pixels), the PDFs obtained 
with the other fields have the same significance. Table B.l 
shows that all the PDFs of Figs. 6 and 7 computed for I — 3 
and 6 pixels are statistically significant as drawn (i.e. with 30 
bins). Those computed for lags I — 12 pixels in the IRAM-30m 
maps are significant only for their smooth characteristics, i.e. 
those visible if the PDFs were built with 8 bins only, as shown 
with the column Z max (8). 



